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We calculate the upper critical field, Hc2{T), due to the orbital pair breaking in disordered 
superconductors without inversion symmetry. Differences from the usual centrosymmetric case are 
highlighted. The linearized gap equations in magnetic field, with the singlet and triplet pairing 
. . . channels mixed by impurity scattering, are solved exactly for a cubic crystal. 
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r^ ! I. INTRODUCTION 
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f— N ' Recently, superconductivity has been discovered in a number of compounds lacking inversion symmetry, such as 

^ ; CePtsSi (Ref. [j, Ulr (Ref. [f), CeRhSis (Ref. S), CelrSis (Ref. i), Li2(Pd3-:rPtx)B (Ref. i), and many others. Much 

of the theoretical work in the field has focussed on searching for the features which are specific to noncentrosymmetric 

systems. These include the magnetoelectric effect,^'^'^ a large residual spin susceptibility and reduced paramagnetic 

limiting)ii2iiiSiiiii^ and various novel nonuniform superconducting statesJ^'i^ii^'i^ 

fj ■ In this paper we study the effects of the absence of inversion symmetry on the upper critical field, Hc2{T), at 

I , arbitrary temperature. We assume the pairing to be of the Bardeen-Cooper-Schrieffer (BCS) type, and include 

O , only the orbital pair breaking. The main qualitative difference from the centrosymmetric case is that the spin- 

^ _ orbit (SO) coupling of electrons with the crystal lattice changes the nature of single-electron states, lifting spin 

"^ ' degeneracy of the energy bands. Then even scalar impurities can mix the singlet and triplet channels in the Cooper 

pair propagator, thus making the theory considerably more complicated. The derivation of the Hc2 equations for 

arbitrary noncentrosymmetric crystal symmetry is presented in Sec. |TT] below, with some of the technical details 

G ' relegated to Appendices El and [B] In Sec. Illli we apply the general equations to a cubic superconductor with the 

1-^ , point group G = O. Assuming that both the band structure and the SO coupling are fully isotropic, we are able to 

r^ • exactly solve the coupled equations for the singlet and triplet channels, obtain the Hc2 equation in a closed form, and 

O \ derive analytical expressions for the upper critical field in the "dirty" limit. This isotropic model clearly shows the 

CJ . deviations from the usual, i.e. centrosymmetric BCS, case, for which the upper critical field was calculated in the 

[ classic papers by Helfand, Werthamer, and Hohenberg in 1960s (Refs. llTlilSh . Sec. IIVI contains a discussion of our 

,—1 . results. 

^ ' The magnetic phase diagram of noncentrosymmetric superconductors has been discussed previously in several works. 

'""' , The upper critical field for a clean three-dimensional Rashba superconductor was calculated in Ref. [l5|, while the 
'^ • effects of disorder in the Ginzburg-Landau regime were studied in Ref. 119!. Two-dimensional case, in which only the 
parama gne tic pair breaking is present, was considered in Ref. 20'- Recently, Hc2 at all temperatures was calculated 
in Ref. l2ll neglecting the impurity-induced triplet channel in the pair propagator in the limit when the SO band 
r^ ' splitting is small compared with the Fermi energy. In this paper, we relax this last condition and include both the 
singlet and triplet channels. 

Throughout the paper we use the units in which h = ks = I- 
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II. DERIVATION OF Hc2 EQUATIONS: GENERAL CASE 
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5t 1 Let us consider a noncentrosymmetric superconductor with the Hamiltonian given hy H = Hq + Himp + Hint ■ The 

first term. 

Ho = ^[eo(fc)(5a/3 + 7(A;)o-c,/3]4a"fc/3, (1) 

k 

describes non-interacting electrons in the crystal lattice potential, where a, (3 =t,i are spin indices, eo(^) is the 
quasiparticle energy counted from tp, and <t are the Pauli matrices. In Eq. ((T|) and everywhere below, summation over 
repeated spin indices is implied, while summation over space and band indices is always shown explicitly. The second 
term in Eq. ([T|), with 7(A;) = — 7(— fc), describes a Rashba-type (or antisymmetric) SO coupling of electrons with 
the crystal lattice. ^^ In addition, there might be a usual (symmetric) SO coupling, present even in centrosymmetric 
crystals. If the latter is included, then a, (3 in Eq. ([T|) should be interpreted as pseudospin projections. Diagonalization 



of Hq yields two non-degenerate bands labelled by the helicity A = ±: 

6(fc) = eo(fc) + A|7(fc)|. (2) 

The Fermi velocities in the two bands are given by v\{k) — d^x/dk. The Fermi-level densities of states are defined 
in the usual way by Nx = "^"^^ X]fc "^KaC^)] C^" is the system volume), and the difference between N^ and iV_ is 
characterized by a parameter 

If the SO coupling is small compared with the Fermi energy, then 5 ^ 0{EsoI^f): where Eso = 2maxfc |7(fe)| is a 
measure of the SO band splitting. 

Scattering of electrons at isotropic scalar impurities is introduced according to 

Hrmp = J d^r U{r)^Pi{r)^P^{r). (4) 

The random potential U{r) has zero mean and is characterized by the correlator {U{r)U{r')) — ni„ipUQS{r — r'), 
where riimp is the impurity concentration and Uq has the meaning of the strength of an individual point-like impurity. 
The field operators are given by V'a(^) = V^^/^ J2k ^^''^'^ka- 

Neglecting the paramagnetic pair breaking, which is a good assumption in many bulk noncentrosymmetric materials, 
the effect of a uniform external magnetic field H is described by the Peierls substitution!^ 

h^eo{K)+j{K)& + U{r), (5) 

where K = — iV + (e/c)A{r), and e is the absolute value of the electron charge. 
We describe the pairing interaction by a BCS-like Hamiltonian: 

H,,,t = -V J d'r^\{r)i:l{r)i:^{r)^^{r), (6) 

where V^ > is the coupling constant. In this model, the superconducting order parameter is represented by a single 
complex function ri{r), see Ref. 1241 . The critical temperature at a given field, or inversely the upper critical field, 
Hc2{T), at a given temperature, is found from the condition that the linearized gap equation 



i-Tj^'lK) 



V 



ri{r) = (7) 



has a nontrivial solution. Here u;„ = (2n -I- l)7rT is the fermionic Matsubara frequency, the prime in the second term 
means that the summation is limited to |w„| < Wc, where lOc is the BCS frequency cutoff, and the operator X(w„) is 
defined by the following kernel: 

X(r, r'; ^„) = i( t^9^G{r, r'; c^„)gG^(r, r'; -c^„))„„p, (8) 

where g = ia2- The angular brackets denote the impurity averaging, and G{r,r']ujn) is the Matsubara Green's 
functions of electrons in the normal state, which satisfies the equation 

(iw„~/i)G(r,r';w„) = <5(r-r'), (9) 

where the single-particle Hamiltonian h is given by expression ([5|). 

At zero field, Eq. ^ yields the following expression for the average Green's function: 

Goik.un) = J2 nA(fe)GA(fc,c^„), (10) 

A=± 

where 

n.(fe) = ^^4^ (11) 
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FIG. 1: Impurity ladder diagrams in the Cooper channel. Lines with arrows correspond to the average Green's functions of 
electrons, g = ia2, and the impurity (dashed) lines are defined in the text, see Eq. (fT5)) . 

are the band projection operators (7 = 7/I7I), and 

Gx{k, ^„) ^ -. \ . , (12) 

are the electron Green's functions in the band representation. Here ^x{k) is the quasiparticle dispersion in the Ath 
band, see Eq. ^, T = l/2r is the elastic scattering rate, r = {2TTnimpUQ Np)~^ is the electron mean free time due to 
impurities, and 

N. = ^^. (1.) 

The impurity average of the product of two Green's functions in Eq. ([5]) can be represented graphically by the 
ladder diagrams, see Fig. [1] We assume the disorder to be sufhciently weak for the diagrams with crossed impurity 
lines to be negligible, see Ref. l25l In order to solve Eq. ([7]) at nonzero field, we introduce an impurity-renormalized 
gap function £)(r,a;„), which a matrix in the spin space satisfying the following integral equation: 

D{r, tj„) = ri{r)g + -UimpU^g / cfr' trg^G{r, r'; w„)i)(r', cj„)G^(r, r'; -cj„) 

+ ^n.mpC/o9 / d^r' tr g^G{r,r';uJn)D{r',uJn)G'^{r,r';-ujn), (14) 

where G{r, r'; LOn) are the disorder-averaged solutions of Eq. ([9]). The above equation can be easily derived from the 
impurity ladder diagrams in Fig. [U by representing each "rung" of the ladder as a sum of spin-singlet and spin-triplet 
terms: 

TI'irnpUQOp^i/Opu — 'i^^^imp'-J Q 9fj,p9au ~^ T^^imp'J ()Q f.ipQau': \^'^) 

where g = ia(j2. 

Seeking solution of Eq. (fTi|l in the form 

b{r, w„) = do(r, w„)g -(- d(r, w„)97 (16) 

we obtain a system of four integral equations for rfa(r, a;„), where a — 0, 1,2,3: 

3 

fc=0 

Here the operators yah{^n) are defined by the kernels 

:y,b(r,r';c^„) = _L- trgtG(r,r';a;„)gfcG^(r,r'; -c^„), (18) 

with go = g, and gi = gi for i = 1,2,3. We see that, in addition to the spin-singlet component do(T',Ci;„), impurity 
scattering can induce also a nonzero spin-triplet component d{r,ujn)- The gap equation ([7]) contains only the singlet 
component: Using Eqs. (fT7|). we obtain: 

NpV "^ ' ^ T ^ ' 

n 

We would like to note that the triplet component does not appear in the centrosymmetric case. Indeed, in the absence 
of the Zeeman interaction the spin structure of the Green's function is trivial: Gapir, r'; ujn) = SapG{r, r'\ a;„). Then 
it follows from Eq. P^ that yab{'-^n) = '^ab3^(w„), therefore do = (1 ~ ^y)^^^] and d = 0. 



The next step is to find the spectrum of the operators yab{i^n)- The orbital eflFect of the magnetic field is described by 
a phase factor in the average electron Green's function: G{r, r'; ujn) = Gq{i — r'; (jjn)e^'^^'^'^ \ where Gq is the average 

Green's function in the normal state at zero field, ip{r, r') — (e/c) / A{r)dr, and the integration is performed along 
a straight line connecting r and r' (Ref. 25). The "phase-only" approximation is legitimate if the temperature is not 
very low, so that the Landau level quantization can be neglected. Using the identity e^"^'^'"''' ^rjir') = e^*^'' '' ^^ri{r)^ 
where D = — iV + (2e/c)A, we obtain: 

yab{'jJn)=yab{q-,UJn)\^^jj, (20) 

where 

yab{q,ujn) ^ ^-j^ -——tTglGo{k + q,uJn)gbG'^{-k,-uJn), (21) 

Substituting here the Green's functions PH)) and calculating the spin traces, we obtain for the singlet-singlet term: 

yoo{q,^n) = 7;Y]Px(] 1 , T^ , . 77T : 777/ : (22) 

where 

P^^^ = 1±S (23) 

are the fractional densities of states in the two bands, and ((...)) a denotes the Fermi-surface averaging in the Ath 
band. Similarly, for the singlet-triplet mixing terms we obtain: 

yo.iq,u;,,)=y,oiq,u;.n) = lj2^pU \ l^r^-^'^M • 7^) ' ^^4) 

2-^ \|cj„| +r + zt;A(fc)gsignw„/2/^ 

We see that the mixing occurs due to the SO coupling and vanishes at 7 -^ 0, when p+ = p_ = 1 and v^ = v_ = Vp. 
Finally, the triplet-triplet terms can be represented as follows: 

y^M,Lo^) =y\f {q,u^) + y\f {q,Lu^), (25) 



where 



and 



3^,, (9,Wn)-^2.PA(r-^— ^T-— 7TT-— — -7^) , (26) 



2^ \ |w„| +r + it;A(fe)q'signtj„/2 



y\f{q.^n) = l^^^^Y^ j -7^{Sij - J^^j " iXeijai)Gx{k + q,u;n)G-x{-k,-uJn). (27) 

The singlet impurity scattering channel, which is described by the first term in expression (|15p . causes only the 
scattering of intraband pairs between the bands. In contrast, the triplet impurity scattering can create also intcrband 

- (2") 

pairs, which are described by y^j ■ It is easy to show that if the SO band splitting exceeds both ujc and F, then the 
second (interband) term in Eq. ([^S]) is smaller than the first (intraband) one, see Ap pendix \X\ Note that in real 
materials, Eso ranges from tens to hundreds meV, see Ref. 26 for CePtsSi, and Ref. 127 for Li2Pd3B and Li2Pt3B. On 
the other hand, there is still considerable uncertainty as to the values of lUc, especially in hcavy-fermion compounds, 
such as CePtaSi. The typical energy of phonons responsible for the pairing in Li2Pd3B was estimated in Ref. 27 to 
be 20 meV, while the SO band splitting is 30 meV (reaching 200 meV in Li2Pt3B). 

The critical temperature of the phase transition into a uniform superconducting state at zero field can be found by 
setting q = in the above expressions. According to Eq. ([2^. the singlet and triplet channels are decoupled. Then 
it follows from Eqs. ((22|) and ([T7| that do{uJn) = (1 -t- F/jWnDTy. Substituting this into Eq. (fTQ]) . we obtain: 



NfV 



-^E'r^-0, (28) 



which yields the superconducting critical temperature: 

r,o = ^o;ee-i/^-^, (29) 

TT 

where C ~ 0.577 is Euler's constant. We see that there is an analog of Anderson's theorem in noncentrosymmetric 
superconductors with a BCS-contact pairing interaction: The zero-field critical temperature is not affected by scalar 
disorder ji^ 

In the presence of magnetic field, neglecting the interband contributions to the triplet pair propagator, we obtain: 

.-, . ^ IV- / AA,a(fc)AA,fc(fc) \ ,_. 

2^ \\uJn\ + T + ivx{k)qsignuJn/2/ ^ 
where 

Next we use in Eq. ([501) the identity x^^ — L due~^"', and make the substitution q ^ D, see Eq. ([^(7)) . in the 
exponent to represent 3^q6(w„) as a differential operator of infinite order: 

yabi^u) = 2 / ^" e-"(l-"l+r) ^ p,Of , (32) 

Jo ^ 

where 

Of - (AA.a(fc)AA,fc(fc)e-'"''^('=)-°^'snc^„/2\ (33) 

In order to solve Eqs. (fT7)) . with the operators !Va6(wn) given by expressions ([5^ . we follow the procedure described 
in Ref. [l3- We choose the z-axis along the external field, so that H = Hz, and introduce the operators 

a±=£H^^^4^' "3-^/fI)., (34) 

where in = y^c/eH is the magnetic length. It is easy to check that a+ = a_ and [a_,a+] = 1, therefore a± have 
the meaning of the raising and lowering operators, while 03 — Og commutes with both of them: [03, a±] = 0. It 
is convenient to expand both the order parameter 77 and the impurity-renormalized gap functions da in the basis of 
Landau levels \N,p), which satisfy 



a+\N,p)^VN+l\N + l,p), a^\N,p) = VN\N-l,p), a3\N,p) = p\N,p), (35) 

where A'^ = 0, 1, ..., and p is a real number. We have 

V{r) = ^w,p(r|Ar,p), da{r,iun) = ^ d^,p(a;„)(r|A^,p). (36) 

N,p N,p 

According to Eqs. (fTT]) . the expansion coefficients satisfy the following algebraic equations: 

y^ ISabSNN'Spp' -T{N,p\yab{(^n)\N',p') d%,p,{uJn) = 6aOriN.,p- (37) 

N',p',b 

Substituting the solutions of these equations into 

1 _ _^y^.<pK)-W.p^p^ ^33^ 



:1lN,p 



NfV ' '' "^ r 

see Eq. (J19p . and setting the determinant of the resulting linear equations for r\N,p to zero, one arrives at an equation 
for the upper critical field. 



III. CUBIC CASE 



In the general case, i.e. for arbitrary crystal symmetry and electronic band structure, the procedure outlined in the 
previous section does not yield an equation for Hc2{T) in a closed form, since all the Landau levels are coupled, and 
one has to diagonalize infinite matrices. In order to make progress, we focus on the case of a noncentrosymmetric cubic 
superconductor with the point group G = O, which describes, for instance, the crystal symmetry of Li2(Pdi_2;,Pta:)3B. 
The simplest expression for the SO coupling compatible with all symmetry requirements has the following form: 



7(A;) = 7ofc, 



(39) 



where 70 is a constant. We assume a parabolic band: eo(fc) = k'^ /2m* — ep, where to* is the effective mass, 
tp = k'^l2m* , and fco is the Fermi wave vector in the absence of the SO coupling. The band dispersion functions are 
given by 



ix{k) = 



/r2 



2to* 



A|7o|fc, 



(40) 



so that the SO band splitting is isotropic and given by Eso 
coupling strength by a dimensionless parameter g — Eso/ 4:6 f 

Xg), the Fermi velocities are the same: vx{k) = vpk, where vp 



«F,A 



fco(Vl^ 



g- 



2|7o|^o- It is convenient to characterize the SO 
While the two Fermi surfaces have different radii: 

fco\/l + g^/m*. For 



the parameter 5, which characterizes the difference between the band densities of states, see Eq. (j3]), we have 

\5\ = 2g^l + g^/{l + 2g^). We assume that 

5c « |<5| < 1, (41) 

where 5c = max(tt'c,r)/ei7' ^ 1. While the first inequality is equivalent to the condition Eso S> max(ti;c,r), which 
ensures the smallness of the interband contribution to the Cooper impurity ladder (see Appendix [X| , the second one 
is always satisfied, with |i5| -^ 1 corresponding to the rather unrealistic limit of extremely strong SO coupling, g — > 00. 
In order to solve the gap equations, we make a change of variables in the triplet component: 



d±{r,ujn) 



di{r,LL!.n)±id2{r,uj„) 
V2 



Then, Eqs. (fT7|) take the following form: 



/ 



1 



-ryos 1 - ry-is 
-ryo+ -r%+ 
V -rj^o- -ry-i- 



-ryo- -ryo 

-ry^^ -rys 

1 - rz -rz_ 

-ri_ 1 - ri / 




(42) 



where 



3^0+ = 



ypi ± iyo2 

V2 



y3±^ 



yi3 ± »y23 

^/2 ■ 



Z = 



3^11 + 3^22 



z± = 



3>11 ± 2Z3>12 - 3^22 



with yab = yba given by Eqs. ^^. 

According to Sec. [Til one has to know the matrix elements of the operators yabi'^n) in the basis of the Landau 
levels \N,p). After some straightforward algebra, see Appendix [Bl we obtain the following expressions for the nonzero 
matrix elements: 

{N,p\yoo{^n)\N,p) ^ y^j^JUn), 
{N,p\%3{lo^)\N,p) = y%^Jto^), 

{N,p\Z{uJn)\N,p) = ZN^p{uJn), 



where 



Jo Jo 

/•OO />1 

Jo Jo 



2Af,p 



2 Jo Jo 



(43) 
(44) 
(45) 
(46) 



V — (vp sigiiUJn/2iH)u, and Ln{x) are the Laguerre polynomials of degree N. Similarly, we obtain: 

{N,p\yo-itOn)\N + l,p) = {N + l,p\yo+{Lu,,)\N,p) = y^^piL,^), 
{N,p\%_{lo^)\N + l,p) = {N + l,p\%+iuo,,)\N,p) = yljun), 
{N,p\Z.{uJn)\N + 2,p) = {N + 2,p\Z+iuj„)\N,p) = SnA^u), 



where 



yN,pi^n) = -iS 



VN-pi^n) 



^2{N+l)Jo 



d^e-n(|^„|+r) / dsv{l~ s')cosipvs)e-''"'^'-''^/'L'^^\v\l- s% 



^2{N + l)Jo 



ZN,p{^n) — — 



-w, 



2^{N + l){N + 2)Jo 



(47) 
(48) 



rf^e-«(KI+r) / ds„s(l_52)giJ^(p^^)g-^-=(l-s=)/2^(l)[„2(^_^2)] 



and LjY (x) are the generalized Laguerre polynomials. 

It follows from the above expressions that the Landau levels are decoupled, and for ri{r) — r]{r\N,p) [rj is a constant) 
the solution of Eqs. (|42|) has the following form: 



do(r,u;„) \ / d%Aun){r\N,p) \ 
- ■ ■ ' <.p(u;„)(r|7V,p) 

rf+_p(a;„)(r|7V + l,p) 
Vd^,p(^„)(r|7V-l,p)/ 

For given N and p, the coefficients are found from the equations 

y^ Mab{N,p;UJn)d%p{uJ„) ^ SaOVi 



d^ir^ujn) 
d+(r,w„) 
d_(r,w„) 



(50) 



(51) 



6=0,3, ± 



where 



M{N,p;uJr, 



~^yN,p 1 ~ -i^^^^Tv,;) ~^yN.p ~^yN-i.p 



(52) 



Substituting the solution of Eq. (|5ip in Eq. (|38p . and using Eqs. (pS)) and (P^)) to eliminate both the frequency cutoff 
and the coupling constant, we obtain an equation implicitly relating the magnetic field and the transition temperature 
at given TV and p: 



1 Zi£= yv^j 1 [M-\N,p;u;^)]oo-l 

^ T ^ 1^\ I,,, I r 



(53) 



The upper critical field, Hc2(T), is obtained by maximizing the solution of this equation with respect to both N and 
P- 



Note that the matrix elements of A4 which are responsible for the singlet-triplet mixing, i.e. y}f„ 



y%. 



and 



y%-ip, are all proportional to (5, see Eqs. dH]) and (|T7)) . Therefore, at \5\ <C 1 the singlet and triplet channels are 

effectively decoupled. Neglecting the corrections of the order of (5^, we obtain from Eq. ([5T|) that [M~^{N,p\u)n)\oQ = 
(1 — Fy™ )~^. Substituting this into Eq. ((53|) . we recover the Helfand-Werthamer expressions,— with the maximum 
critical field corresponding to A^ = p = at all temperatures. Thus, in the weak SO coupling limit the absence of 
inversion symmetry does not bring about any new features in Hr^lT), compared with the centrosymmetric case (as 
long as the paramagnetic pair breaking is not included, see Ref. l2ll ). 

A. "Dirty" limit at iV = 0, p = 

At arbitrary magnitude of the SO band splitting, the singlet-triplet mixing makes the Hc2 equation in noncen- 
trosymmetric superconductors considerably more cumbersome than in the Helfand-Werthamer problem, even in 
our "minimal" isotropic model. It is even possible that, at some values of the parameters, the maximum crit- 
ical field is achieved for A^ > and p ^ Q, the latter corresponding to a disorder-induced modulation of the 
order parameter along the applied field. Leaving investigation of these exotic possibilities to future work, here 



we just consider the case N = p 

[^-i(0,0;a;„)]oo = (l-rzi,o)/[(l- 
ature, magnetic field, and disorder: 



0. Then it follows from Eqs. ^ and §1 



that (Iqq — 



^0,0 



0, and 



r2/o'o)(l — Ezi^o) — r^(yo q)^]. It is convenient to introduce the reduced temper- 



t = 



T 

Trn 



2H 



-TT, C 



ttT, 



cO 



where Hq = ^q/tt^^q, $0 = nc/e is the magnetic flux quantum, and ^0 = vf/^ttTco is the superconducting coherence 
length. In these notations, Eq. ((53|) yields the following equation for the upper critical field hc2{t): 



In- 



1 



n>0 



2„2 



2n + l il~Cw^){l-CPn) + C^S^ql 



(54) 



where 



Jo Jo 

Pnit,h)^ / dpe-'^-P / ds^— 
Jo Jo ^ 



2' ^ 



-hp^{l~s^)/4 



(55) 



where a)„ — {2n + l)t + (. 

In the clean limit, i.e. at C — > 0, or if the SO band splitting is negligibly small, i.e. at (5 ^ 0, one recovers from Eq. 
(I54p the Helfand-Werthamer equation for a centrosymmetric superconductor. Thus the absence of inversion symmetry 
affects the upper critical field only if disorder is present. One can expect that the effect will be most pronounced in 
the "dirty" hmit, C 3> 1. [Note that, according to Eq. (|4T|l . the disorder strength should satisfy C ^ {^f/Tco)\S\-] We 
shall see that in this limit hc2 scales as C, which allows one to use the Taylor expansions of the exponentials in Eqs. 
dST 



Wnit,h) 



I 



I- 



3(1), 



7,2 



Pn{t,h) 



1 

3a)„ 



1- 



6/i 



5a); 



qn{t,h) 






Using the fact that the main contribution to the Matsubara sum in Eq. ([M|) comes from {2n + l)t ^ (, we arrive at 
a well-known universal equation, which describes the magnetic pair breaking in superconductorsi^ 



where "^(x) is the digamma function, and 



t \2 t ' 



2 + S^ 



12c 



(56) 



(57) 



characterizes the pair-breaker strength. Note that the corresponding expression in the centrosymmetric case is differ- 
ent: acs = h/6C (Ref. 1171 ). Analytical expressions for the upper critical field can be obtained in the weak- field limit 
near the critical temperature: 

and also at low temperatures: 

hc2\t=o = I^C- (59) 

We see that the SO band splitting in the noncentrosymmetric case enhances the orbital pair breaking. 

IV. CONCLUSIONS 

We have derived equations for the upper critical field in noncentrosymmetric superconductors, assuming a BCS- 
contact pairing interaction and orbital pair breaking. In a cubic crystal (the point group G — O), in which both 
the electron dispersion and the SO coupling are isotropic, the gap equations are shown to be diagonal in the Landau 
level basis, with the singlet and triplet channels in the pair propagator mixed together. For the order parameter 
corresponding to the lowest Landau level, without any modulation along the applied field, we have derived the -ffc2 
equation in a closed form and solved it in the "dirty" limit, in which the effects of the absence of inversion symmetry 
are expected to be most pronounced. The effect on the upper critical field of the singlet-triplet mixing, which is 
responsible for the deviations from the Helfand-Werthamer theory, is found to be proportional to S^. 

Application of our theory to real noncentrosymmetric superconductors of cubic symmetry, such as Li2(Pd3_2;Pt2,)B, 
is complicated by the fact that the Fermi surfaces as well as the SO band splitting are strongly anisotropic. Using 
the maximum values of the SO band splitting from Ref. [23, one can estimate the corrections to Hc2 (T) due to the 
singlet-triplet mixing to be of the order of several percent. 
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APPENDIX A: INTERBAND VS INTRABAND CONTRIBUTIONS 

In this appendix, we estimate the relative magnitudes of the intraband and interband contributions to the triplet 
pair propagator, Eq. (j25p . in the limit when the SO coupling is strong compared with both the cutoff energy ujc and 
the elastic scattering rate F. Let us consider an isotropic band with -f{k) — 70A; in a cubic crystal. Neglecting for 
simplicity the differences between the densities of states and the Fermi velocities in the two bands: p+ ^ p- = I and 
V-^- = v^ = vp, and setting q = 0, we obtain from Eqs. ((26l) and ([27]) : 

-r,(l)/ n \ (7J7i)fc ^ij ,, f sr 



and 



2 4" \ |w„| +T + zA|7| sign w„ /^ 3(|c^„| + F)(l 



A 



yij {q - o,a;„) - ^ E ( rTjTTrr-rnrf— — - ).. - qm^.j 4-nn -j-r-2^ = ^™*-K)'5y' 



where r(a;„) = Eso/'^i\^n\ + T). Due to the BCS cutoff, the maximum value of w„ is equal to lUc, therefore r„ 
i?So/max(u;c,F) :g> 1. From this it follows that 



y inter y-^n) 

max — - = - 



2 

rain 



max(cjc, F^"' 



Eso 



<1. 



Therefore the interband contribution is small compared with the intraband one, at all Matsubara frequencies. 
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APPENDIX B: CALCULATION OF {N,p\yab{uJn)\N' ,p') 

The operators yab{i^n) are given by expressions ([32)) . For a spherical Fermi surface and 7(fc) — 7ofc, we obtain 
froniEq. §^: 

of = i / d6l sin 6*6-'™='™"^ / ^$f (6»,,^)e-™('=""*''++'="*°")'''"^ (Bl) 

2 Jo Jo 27r 

where v = (u_f signu;„/2^/f)M, and ^f'{9,(j)) = AA,a(fc)AA,b(fc), with AA,o(fe) = 1 and AA,i(A;) = Xki for i — 1,2,3, see 
Eq. ([5T|) . Using the well-known operator identity e"^"*"^ = e~["^''^l/^e'^e^, which holds if the commutator of A and B 
is a c- number, and expanding the exponentials in powers of a±, we obtain: 

yab{oJn) = -y2p>^ dwe-"(l'""l+^M d0sin0e-™'^^^°''«e-("'/2)sm^e^ab^0)^ (g2) 

4 ^ Jo Jo 

where 






n,m— 



27r 



(B3) 



Below we perform the detailed calculations for Jpoo and 3^o- = (iVoi ^ iyo2)/V^- Other matrix elements can be 
considered in a similar fashion. 

ypp : Since ^^{6, </>) = 1, the (/)-integral on the right-hand side of Eq. (jB3p is equal to (5„m, and 

^A (^) - 2^ r-jT^ a+a-- 

»i=o ^ ■' 

It is easy to show, using Eqs. p5|) . that a"a" |Af,p) = [Nl/{N — n)\]\N,p) for n < N, and zero otherwise. Therefore, 

^ ATI 

where Lm{x) is the Laguerre polynomial of degree N . Substituting this into Eq. (jB2p . using the fact that p+ +p_ = 2, 
and introducing s = cos0, we obtain: {N,p\yoo{i-jJn)\N,p) = j/™p(w„), where y^p{ujn) is given by Eq. (|43|l . 
j)o-: Since «f>A"(0, 0) = Asin6'e~*'^/\/2, one has m = n + 1 on the right-hand side of Eq. ([BSJ) . and 



^r(^) = A^E 



smf v-^ —iwsmi 



\2n+l 



n ^n+1 



V2 ^ ?i!(n-M)! 

^ n— 



-a I a 



V^ n=0 ■-■V < 

Using a!j:a"+^|iV +l,p) = ^/N + 1[N\/{N - n)!]|iV,p) (if n < N, zero otherwise), we obtain: 



t2/) ^ AT\ ('_,,2 oiTi2 fl^n 



^ n— / \ / 



, ^ 7,sin2flr.W^i^c;r.2, 

v/2(]vTT) "^ 



-a - 7,sin2fl/-.(.V(„2sin2fl)|A^ J,) 



where L^ {x) are the generalized Laguerre polynomials: 



TV 



i^"^(-) = E 



(iV + a)! (-a;)" 



(iV-n)!n!(n-f a)!' 

see Ref. [2^ (the Laguerre polynomials are recovered by setting a ^ 0). Substituting this into Eq. (|B2p and using 
p+- p- ^ 2S, we obtain: (iV,p|J^o_(w„)|Ar + l,p) = y^p(w„), where y%^p{uJn) is given by Eq. (gT]). 
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